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Abstract 

We study the thermodynamic properties of an ideal gas of fermions in a harmonic 
oscillator confining potential. The analogy between this problem and the de Haas- 
van Alphen effect is discussed and used to obtain analytical results for the chemical 
potential and specific heat in the case of both isotropie and anisotropie potentials. 
Step-like behaviour in the chemical potential, first noted in numerical studies, is 
obtained analytically and shown to result in an oscillatory behaviour of the specific 
heat when the particle number is varied. The origin of these oscillations is that 
part of the thermodynamic potential responsible for the de Haas-van Alphen-type 
effect. At low temperatures we show analytically that there are significant deviations 
in the specific heat from the expected linear temperature dependence, again as 
a consequence of the de Haas-van Alphen part of the thermodynamic potential. 
Results are given for one, two, and three spatial dimensions. In the anisotropie case 
we show how the specific heat jumps as the ratio of oscillator frequencies varies. 
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1 Introduction 

The experimental realization [1] of Bose-Einstein condensation in confined 
gases of atoms at low temperatures has been the stimulus for a wide range 
of experimental and theoretical investigations. (See Ref. [2] for a review.) By 
ingenious experimental techniques it is possible to prepare the atoms as either 
bosons or fermions. There have now been a number of experiments on trapped 
gases of fermions, as well as mixtures of bosons and fermions [3]. 

For a gas of fermions, it is a good approximation to model the system as 
a collection of non-interacting particles obeying Fermi-Dirac statistics con- 
fined by a simple harmonic oscillator potential. Unlike the case of bosons, the 
dominant s-wave scattering channel is suppressed making the effects of inter- 
actions less important [4]. With this simple theoretical model for the system, 
it is easy to obtain exact and simple expressions for the single-particle energy 
levels using the standard quantum mechanical result for the simple harmonic 
oscillator. This can serve as a starting point for a theoretical analysis of the 
thermodynamic behaviour of the system. 

Because the energy levels for a set of particles (bosons or fermions) confined 
by a simple harmonic oscillator potential form a discrete set, in order to calcu- 
late the thermodynamic properties we must be able to perform sums over the 
quantum numbers labelling the states; this can lead to difliculties in evalua- 
tion. By restricting the parameters of the system (temperature and oscillator 
frequencies for example) it may be possible to argue that the sums can be 
approximated with integrals, an approximation that renders the computation 
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of analytical results easier. We will refer to the approximation of sums with in- 
tegrate as the continuum approximation, since it is equivalent to regarding the 
discrete energy spectrum as continuous. This approximation was first applied 
to the Fermi gas in the now classic paper of Butts and Rokhsar [4] . Even when 
one goes beyond this approximation, as we will show below, the continuum 
approximation can stili be very accurate in some temperature regimes, and 
can capture the leading order behaviour if the temperature is not too low. 

The first study that showed there might be more features present beyond 
what is obtained by the continuum approximation was given by Schneider 
and Wallis [5]. These authors computed the thermodynamic expressions nu- 
merically, and their results showed that there were a number of step-like fea- 
tures present, not predicted by the continuum approximation. In particular, 
the chemical potential, which is a smooth function in the continuum approxi- 
mation, has a series of step-like jumps when a numerical analysis of the exact 
result is performed. This in turn can lead to step-like, or oscillatory, behaviour 
of other thermodynamic properties, such as the specific heat. Another impor- 
tant consequence of the work presented in Ref. [5] is that at sufnciently low 
temperatures there are significant deviations in the specific heat from the lin- 
ear temperature dependence predicted [4] by the continuum approximation. 
The main purpose of the present paper is to analyze the situation analytically 
to obtain the step-like behaviour of the chemical potential, and to apply this 
method to the evaluation of the specific heat. 

In addition to the two studies already mentioned [4,5], there have been a 
number of other related theoretical works. Ref. [6] generalized the calculations 
of Ref. [5] to a gas with two spin states with an interaction between the two 
states. (We will make some brief comments on the effects of interactions in 
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Sec. 6.) A number of authors [7,8,9,10,11,12,13,14] have obtained analytical 
approximations in a number of cases using path integrai, Green's function, 
or density matrix methods. The approach described in the present paper was 
used to study partially confined gases in Ref. [15] 

In the present paper we will concentrate on the behaviour of the chemical 
potential and sperine heat, and will use the evaluation of the thermodynamic 
potential as the basic starting point. The key observation is that for a simple 
harmonic oscillator confining potential, the evaluation of the thermodynamic 
potential is mathematically very similar to that of a gas of charged fermions in 
a Constant magnetic field. Landau [16] was the first to show that the thermody- 
namic properties of electrons in a Constant magnetic field undergo oscillations 
whose period is determined by the inverse of the magnetic field strength. This 
gave a theoretical basis for the observations of de Haas and van Alphen, now 
referred to as the de Haas-van Alphen effect. (See Ref. [17] for a good dis- 
cussion of the de Haas-van Alphen effect and the important role it plays in 
condensed matter physics.) What we will show here is that for the trapped 
Fermi gas the step-like behaviour found in Ref. [5] is completely analogous 
to the de Haas-van Alphen effect. Just as in the de Haas-van Alphen effect, 
where use of the continuum approximation would miss the oscillations, here 
too we must go beyond the continuum approximation. 

The outline of our paper is as follows. In Sec. 2 we describe the general method 
following the classic analysis of Sondheimer and Wilson [18]. There are some 
differences with the de Haas-van Alphen case considered in Ref. [18] that we 
will describe. We will show how the continuum approximation comes about in 
the method, and even find the next order correction to it, beyond the leading 
order result given in Ref. [4]. In Sec. 3 we analyze the de Haas-van Alphen 
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contribution to the thermodynamic potential for a 3-dimensional isotropie 
harmonic oscillator potential, and obtain approximate analytical results for 
the chemical potential. This is then used to examine the specific heat. Results 
for the 1- and 2-dimensional gases are given in Sec. 4 and the results compared 
with the 3-dimensional case. In Sec. 5 we analyze the situation when the 
trapping potential is anisotropie. Finally, Sec. 6 contains a brief summary and 
a short discussion of the main results. 

2 General method 

We will begin with the grand canonical ensemble and the thermodynamic 
potential Q for a system of fermions with energy levels E n defined by 



We use the usuai notation (3 = 1/T, with T the temperature, and choose units 
such that the Boltzmann Constant is equal to one. Take the Laplace transform 
of f(E), and cali it </?(/?), so that 





n 



where 




(2.2) 



OC' 




(2.3) 



o 



The inverse Laplace transform of Eq. (2.3) reads 




c+ioo 



(2.4) 
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where c is an arbitrary Constant chosen so that the integration path lies to the 
right of any singularities of ip(P). We can now use Eq. (2.4) in Eq. (2.1). The 
sum over the energy levels that results can be related to the partition function 
for the /i — system, defined by 

^) = E^ E "- (2-5) 

n 

We can then write Q defined by Eq. (2.1) as 

c+ìoo 

il = -— J d(3Z{-(3)y{(3) . (2.6) 



C— lOO 



It is necessary here to regard Z(f3) as a function of (5 with f3 viewed as a 
complex variable. Z{—j3) is the result obtained from this complex function by 
analytic continuation of the definition Eq. (2.5) to 9?(/3) < 0. 

We now define Ò(E) to be the Laplace transform of (3~ 2 Z((3). (The factor of 
j3~ 2 is for later convenience [18].) This results in the definition 

oo 

Z((3) = (3 2 JdEe-P E 3(E). (2.7) 
o 

The inverse Laplace transform gives 

e+ioo 

3(£) = — / d(3eP E (3- 2 Z{(3). (2.8) 

c— ioo 

Substitution of Eq. (2.7) into Eq. (2.6) followed by the use of Eq. (2.4) results 



m 



n = -TjdE3(E)^±. (2.9) 
o 



The Fermi-Dirac distribution function F(E) is 
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F(E) = [e^ E -rì + l] 1 (2.10) 
= -T^f(E), (2.11) 



if we use Eq. (2.2). This leads to 

= JdEÒ(E)—F(E), (2.12) 



giving the key starting point in the Sondheimer- Wilson [18] analysis of the 
de Haas-van Alphen effect. The method rests on an evaluation of 3(E) defined 
by Eq. (2.7) or Eq. (2.8) and its use in Eq. (2.12). 

To evaluate "ò(E) in Eq. (2.8) we require the partition function, or at least 
some information about the singularities of (3~ 2 Z(f3) in the complex /3-plane. 
For the case of a simple harmonic oscillator potential this information is very 
easy to obtain. We will consider the D-dimensional oscillator potential 

\/(x) = ìmf:^x J 2 , (2.13) 

with x = (xi, . . . , x D ) the D spatial coordinates. The energy levels are simply 
given by (in h — 1 units) 

D , 1 



J'=l 



where n = (ni, ... , no) and rij = 0, 1, 2, . . . for j = \,...,D. We then have 
Z(f3) in Eq. (2.5) expressed as a product of geometrie series that are easily 
summed to give 

D -pWj/2 

m = n jr ^-y (2.15) 
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Because of the neglect of interactions, and the simple forni of the potential, it 
has been possible to evaluate the partition function in closed forni. 

To obtain ~ò{E) we need to know the singularity structure of f3~ 2 Z(f3). It is 
now obvious from Eq. (2.15) that (3~ 2 Z(f3) is a meromorphic function. There 
is a pole at (5 — of order D + 2 as well as a series of poles along the 
positive and negative imaginary /3-axis at (3 — (5^ = 2'irikj/ujj with kj = 

±1,±2, The order of these poles depends on the relative ratios of the 

harmonic oscillator frequencies. If none of the oscillator frequencies are rational 
multiples of any of the others, ali of the poles, apart from the one at (5 — 
0, will be simple. However, if some of the oscillator frequencies are rational 
multiples of the others, some of the poles away from the origin will be higher 
order. This makes a determination of the residues of poles along the positive 
and negative imaginary /3-axis for a general harmonic oscillator potential an 
inelegant treatment of special cases. We will return to this problem later. 

We can write "ò{E) in Eq. (2.8) as 



where we dose the contour in the left hand side of the complex piane and use 
3o(E) to denote the contribution to the integrai from the pole at (5 — 0, and 
~òr{E) the contribution coming from the rest of the poles along the positive 
and negative imaginary /3-axis. When Eq. (2.16) is used in Eq. (2.12) we obtain 
the thermodynamic potential as 



3(E)=ò (E)+òr(E), 



(2.16) 



f2 — SÌq -|- fì r 



(2.17) 
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in an obvious way. As in the de Haas-van Alphen effect, the oscillations in 
thermodynamic quantities will come from Q r , but we will first examine Qq. Q r 
will be evaluated in Sec. 3 for the isotropie potential, and in Sec. 5.1 for the 
anisotopic potential. 

To obtain 3o(-^) w e need the residue of [3~ 2 e (5E Z([3) at (3 — 0. For general 
values of D and arbitrary frequencies ujj this is difficult to write down in any 
simple way. We will concentrate on just the cases D — 1,2,3 here, although 
there is no intrinsic reason why the method cannot be extended to other values. 
A straightforward calculation yields 



3o(£) = £-^ forD = l, (2.18) 
3o(£0 = for£ = 2, (2.19) 

ME) = gygoJ I 240 ^ 4 - 120£ V + u\ + ul) + 7(uf + uj\ + cu 3 4 ) 

+10(w?w| + uolul + colui)] for D = 3. (2.20) 



These results can now be used in Eq. (2.12) to obtain 

H = / dE-Ò^E)-F(E) . (2.21) 



As T — * 0, meaning that — > oo, the Fermi-Dirac distribution function 
approaches a step-function, F(E') — > 9 (fi — E), so that -J^F(E) — > —8(E — fi). 
This crude approximation results in 



fio - -3o(/i) , (2.22) 

and is really the first terni in a systematic expansion due to Sommerfeld. (See 
Refs. [17,19,20] for three derivations of the Sommerfeld expansion.) In our case 
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we find 



fio - -3o(/i) 



(2.23) 



Although the general Sommerfeld expansion contains higher order terms, these 
vanish here because 3o(E) is no more than quartic in E due to our restriction 
that D < 3. For arbitrary D, the expansion Eq. (2.23) will contain more 
terms with increasing powers of T. For low values of T these should be less 
important than those indicated in Eq. (2.23) in any case. Using Eqs. (2.18- 
2.20) in Eq. (2.23) results in 



fio: 



+ 



uj'Ì + uj' z 2 - 4vr 2 T 2 ) , for D = 2 , 



6u;ic<;2 24cuia;2 



+ 



-K + ^ + cj 2 -4tt ì T ì ) 



" 576oJicj 2 cj 3 I 7 ^^ + ^2 + <4) + 10(^2 + ^3 + ^M) 
-m-K 2 T 2 (ul + u 2 2 + u)\) + 1127T 4 T 4 Ì , for £> = 3 . 



(2.24) 
(2.25) 



(2.26) 



The thermodynamic properties of the system follow from a knowledge of the 
thermodynamic potential. In particular, the average particle number N is 
given by 



—(E 



(2.27) 



and the internai energy U is given by 



(2.28) 
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In the most physically interesting case we take D — 3. (The cases of D — 1, 2 
will be given later in Sec. 4.) If we temporarily ignore the contribution Q r to 
fi, we have N ~ N where 



,3 



(2.29) 



A* A* /, .2 , , .2 , , .2 ^_2rp2 



6o;ia;2a;3 2Au\U2Uz 



uj{ + lu 2 + lu 2 - An 2 T 2 ) . (2.30) 



For large values of fi, if we keep only the terni in Eq. (2.30) of order fi 3 , we 
find (with N ~ N ), 

/i 3 ~ Qlu^^N , (2.31) 

showing that fi oc N 1 ^ 3 . For large values of iV we will have fi much larger 
than the average oscillator frequency, consistent with the assumptions made 
in deriving these results. It is easy to obtain the next order correction to 
Eq. (2.31), 

,3 r + (4 + 4 + 4-^) 1 (2 . 32) 



[ 4(6A^cjic<J2^3) 2 / 3 

showing that the correction to Eq. (2.31) becomes increasingly unimportant 
for large values of N. If we use Eq. (2.32) and rewrite the result in terms of 
the Fermi energy defined by 

E F = fi(T = 0) , (2.33) 
it is easy to see that 



/^Ml-^|- (2-34) 
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This agrees with the result in Ref. [4] who used the continuum approxima- 
tion and only the T-dependent part of the second terni in Eq. (2.30). We have 
demonstrated here that Eq. (2.34) holds true even when the next order approx- 
imation for the density of states is included in the continuum approximation. 
However, as we will see in the next section there are important corrections to 
this from terms that come from beyond the use of the continuum approxima- 
tion and support the numerical investigations of Ref. [5]. 

We can also calculate the contribution to the internai energy, that we cali 
Uq, coming from Qq. Using Eq. (2.26) in Eq. (2.28) and then eliminating \x in 
favour of iV with Eq. (2.32) results in the following approximation : 

^s^ M+ y ' <2 - 35) 

assuming that N is large. (The intermediate result for Uq in terms of ji can 
be found in Eq. (3.25) below for the isotropie potential.) 

The sperine heat can be found using 

(2.36) 



dT 



N,uj 



The contribution coming from the approximate result in Eq. (2.35) is then 
seen to be 

C ~ ^{Qu^u^^N^T . (2.37) 



As before this agrees with the result quoted in Ref. [4] demonstrating a linear 
dependence on the temperature. Again we will find that at low temperatures 
the part of Q that has not yet been considered, namely Q r , alters this ex- 
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pected behaviour, confirming the numerical results of Ref. [5]. We consider 
the evaluation of iì r and its effect in the next section. 



3 Isotropie 3-dimensional potential 

3. 1 Thermodynamic potential 

We will examine the case of the isotropie harmonic oscillator potential with 

uji = uj 2 = UJ3 = u , (3.1) 

in which case the results of the previous section simplify. Our aim here is to 
calculate Q r that has been neglected up to now from the analysis. To do this we 
need to evaluate the contribution to il coming from the poles of P~ 2 e f3E Z(P) 
along the imaginary /3-axis away from (3 — 0. From Eq. (2.15), because 

-3/3u>/2 

m = (1 - e ^)3 ' ^ 3 - 2 ) 

there will be poles of order 3 at (3 = (3 k defined by 

Pk — — k (3.3) 
u 

for k = ±1, ±2, . . .. It is a straightforward matter to evaluate the residues of 
P~ 2 e P E Z(p) at P — p k and show that their contribution to Eq. (2.8) gives 

^ /rnx ~(-l)W 2l2 4n 2 k 2 E 2 \ , 7rn/x 
UE) = g ^ (6 + n 2 k 2 _ j cos(2^/c) 

+ E [ -^^kE/u). (3.4) 

fc=i 
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It now remains to use this in Eq. (2.12) and to try to extract something useful 
from the result. 



We have 



64tt 4 ^ /e 4 



8tt 3 ^ A; 3 



(3.5) 



where 



dE 6 + 7r z r - 



o- 2 4vr 2 fc 2 £ 2 \ cos(2nkE/uj) 



ex. 

/ 



cu 2 ; cosh 2 [±/3(£-/i)] ' 
sin(27rA;i?/cj) 



cosh 2 ^-^)] 



In these last two integrals we can make the change of variable 



(3.6) 
(3.7) 



(3.8) 



defining a new integration variable 9. The lower limits on the integrals in 
Eqs. (3.6) and (3.7) then become —f3fi/2. If we look at low enough temper- 
atures, specifically T << fi as we have already assumed, then to a good ap- 
proximation (up to exponentially small terms) we can replace the lower limits 
on the integrals defining A k and B k with — oo. The approximate results then 
become 



2 7 de 



- / 



P J cosh 2 # 



2 , 2 1;t 2 /.- 2 // 2 Ki- 2 /.- 2 //^ 



167T 2 fc 2 g 2 

/5 2 cu 2 



6 + vr 2 r - 



/ 



/2ttA;u 4ttà;#\ 
x cos + — — , 

\ LO /JUJ J 

(3 J cosh^r+7 



2 7 de 



20 \ (2ixkfi 4irke\ 
fi + — ] sin I h 



a; 



(3.9) 



(3.10) 
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The integrals in Eqs. (3.9) and (3.10) may be evaluated exactly using residues. 
Ali of the expressions required may be related to the basic integrai 



J co S (ae + b) Tra cos b 

J cosh 2 fl sinh(fa) ' v ; 



— oo 



and derivatives with respect to the parameters a and b. After some straight- 
forward calculation, using the results for Ak and found from Eqs. (3.9) and 
(3.10) as described, it can be shown that 



E 



(-i) fc 



87r 2 /9 3 cc; 2 ^ k 3 sinh 



_ {[8 ^w(?0) 



(3w 



+47r 2 A;^cothf^)+2/3 : 
V Bili ' 



[3u 



2 u 2 



+7i 2 k 2 (4n 2 + f3 2 to 2 - 4/?V)J cos(27rfc/i/cj) 
+4irk/3n (3u + 2rr 2 A;coth (— — ) sm(27rkfx/u) \ , 



(3.12) 



The presence of the trigonometrie functions in Eq. (3.12) is responsible for the 
oscillations that occur in thermodynamic quantities if we go beyond the lead- 
ing order continuum approximation. The presence of the hyperbolic functions 
in Eq. (3.12) with argument 2j 0- means that unless [3uj > 1 these oscillations 
will be suppressed. The oscillatory behaviour should therefore show up for 
T < u and become more prominent as T is reduced. 

The inclusion of Vt r in the expression used for the thermodynamic potential 
will lead to corrections to the thermodynamic behaviour of the system beyond 
what is found using the continuum approximation of Sec. 2. We will look first 
at how the chemical potential is affected. 
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3.2 Chemical potential 



We can define (from Eq. (2.27) 



N r = - ' "'' 



and then use Eq. (3.12) to obtain 



(3.13) 

7> 



^p ò u A ^ smh (^) L L pcj J 

x sm(2nkfi/uj) + 8tc/3/j, coth (— — ) cos(2irk (x/lj) \ . (3.14) 



This will make an oscillatory correction to the contribution N for the average 
number of particles found using the continuum approximation in Sec. 2. 

The chemical potential may be found by solving 

N = N + N r (3.15) 



for fi, where N is given by (see Eq. (2.30) with cu 1 = u 2 = uj 3 = u) 



Obviously the complicated dependence on fi in N r as given in Eq. (3.14) 
renders an analytical evaluation of fi difficult. We can simplify by using the 
assumption that fi is large ( since we have already assumed that (5 fi >> 1 and 
fi^ cu), and keep only the leading terni in fi. This gives 

y „ V^ (-l) fc sin(27rA;^) 
" **k sinh(^) • (3 - 17) 
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The continuum approximation in Eq. (3.16) stili gives the leading contribution 
to N for large fi; however, N r in Eq. (3.17) can be more important than 
the sub-leading term (proportional to fi) in Eq. (3.16) if the temperature is 
low enough. For (3u of the order of 1 or less, the sum in Eq. (3.17) is well 
approximated by simply the first term and it is clear that the oscillations, 
although present, will have a very small amplitude. We would therefore expect 
that for the temperature range T > uj, Eq. (2.31) or Eq. (2.32) would provide 
a good approximation for the chemical potential. However as the temperature 
is reduced, so that (5uj ^> 1, the amplitude of the oscillations coming from 
Eq. (3.17) become more pronounced and must be taken into effect. 

It is possible to find an asymptotic expansion for N r in Eq. (3.17) valid for 
j3uo 3> 1 (or T <C uo). After some calculation, it can be shown that 



with [x] denoting the largest integer whose value is less than or equal to x. 
(Thus < < 1 can be assumed in Eq. (3.18).) Similar asymptotic expansions 
can be found for the sub-leading terms in N r given in Eq. (3.14); however we 
will not give them here. It can be shown by numerically evaluating the sum 
in Eq. (3.17) and comparison with the analytical approximation in Eq. (3.18) 
that Eq. (3.18) does give an accurate result for large values of (5uj. 



N r ~ 




(3.18) 



where 




(3.19) 
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If we take the limit fico — * oo in Eq. (3.18), the result simplifies further to give 

(3.20) 



2co< 



2 1 



co + 2 



i£ _|_ ì 

w ~ 2 



is not 



Strictly speaking, this further approximation is only valid if 
equal to an integer. In this special case, N r in Eq. (3.17) or Eq. (3.18) can be 
seen to vanish, and we must look at the sub-leading contributions to N r that 
follow from Eq. (3.14). 

If we use Eq. (3.20) for N r in Eq. (3.15) along with Eq. (3.16) for N , it is 
possible to show that the solution for \i is given by the simple expression 

- ^ [(67V) 1 / 3 ] + \ . (3.21) 
co z 

The presence of the greatest integer function in Eq. (3.21) leads to the step- 
like behaviour first found in the numerical studies of [5]. Our results provide 
a confirmation of this behaviour by analytical means. For large values of N, 
our result in Eq. (3.21) shows that these steps will occur roughly for iV ~ £ 3 /Q 
where i is an integer. This agrees with the "magic numbers" found by [5] that 
occurred for N = £(£ + 1)(£ + 2)/6 if we make N, and hence t, large enough. 
We have therefore seen how the step-like behaviour of the chemical potential 
comes about in an analytical way, and traced its origin back to the same type 
of terms that are responsible for the de Haas-van Alphen effect. 

In Fig. 1 we show the result for fi/cu plotted as a function of N over a range 
of N. We have taken N large, but the same type of behaviour can be found 
for smaller values as well. The continuum approximation for fi is shown as the 
smooth solid, almost straight, line, and the simple analytical approximation 
in Eq. (3.21) is superimposed on it. The steps occur at the magic numbers as 
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8000 8200 8400 8600 8800 9000 9200 9400 9600 9800 10000 

N 

Fig. 1. (color online) This plot shows /i/uj plotted over a range of N for three 
sample temperatures. The solid, almost straight, line shows the result found using 
the continuum approximation for the particle number resulting in fi/ co = (6iV) 1//3 . 
The step-function superimposed shows the result found from the very low temper- 
ature approximation of Eq. (3.21). The smooth sinuous curves give the results for 
T = O.lco, T = 0.05a> and T = 0.01u>. The last two curves have been displaced from 
the first one for clarity. 



predicted. We have displaced the curves for different values of (5uj away from 
each other for clarity to show the trend towards the step-like behaviour as the 
temperature is reduced. As the temperature is increased, the amplitudes of 
the oscillations decreases as mentioned above. 
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3.3 Specific heat 



We will first calculate that part of the energy that arises from £l r . Define 



U r =-^(f3Q r ) 



(3.22) 



and use Eq. (3.12) for Q r . After some calculation it can be shown that 



7T/1 



3 oo 



E- 



sin 



^3^sinh(20)^V " ) 



37TV ^ (-l) fc COsh(^) / 2 ^ \ 

V^àì si nh 2 (^0) C °H - J 



(3.23) 



where we only include the two leading order terms in \i since it is these two 
terms we will need to calculate the leading contribution to the specific heat. 

The specific heat was defined in Eq. (2.36). It is straightforward to show that 
this expression is equivalent to 

(dU\ fdN\ 
T,u> \9TJ 







(f) 


1> 



(3.24) 



which is more useful for explicitly calculating C . The presence of the second 
term on the right hand side complicates the evaluation of C, but we will obtain 
an expansion for C in powers of \i. 

If we look at the first term on the right hand side of Eq. (3.24) and use 
U = U + U r with U r given by Eq. (3.23) and 



fi 4 ^ 



fìpu,U> 



8lu 3 IQlu 3 



[cu 2 - Att 2 T 2 ) 
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1 



(17u 4 + 40vrVT 2 - 112vr 4 T 4 ) 



(3.25) 



192(W 3 



we obtain 



( 



dU_ 
df 



ì 



ttVT 
2uj 3 



+ 



( 



<9£/ r 



<9T 



) 



(3.26) 



From Eq. (3.23), counting powers of fi, it can be observed that the second 
term on the right hand side of Eq. (3.26) will contain expressions in fi 3 , fi 2 , . . .; 
thus, it appears as if the leading order behaviour of the specific heat will be fi 3 , 
rather than fi 2 as predicted by the continuum approximation. However, we will 
show that the fi 3 part of Eq. (3.26) cancels with a similar expression coming 
from the second term in Eq. (3.24) leaving the overall leading behaviour of 
the specific heat as fi 2 rather than fi 3 . In any case, we must work consistently 
to order fi 2 , and to this order only the two contributions to U r that we have 
written down in Eq. (3.23) are necessary. We will drop ali terms that result in 
explicit factors of fi in C that are of order fi and lower. 

Using Eq. (3.25) we have 



and it is noted that the second term on the right hand side contains terms 
involving fi 3 , fi 2 , From Eq. (3.15) and Eq. (3.16) we have 



with the second term on the right hand side involving fi 2 , fi, . . ., as well as 




(3.27) 




(3.28) 




(3.29) 
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with the second term on the right hand side involving fi 2 , fi,.... A simple 
counting of powers of fi, using the expressions for U r and N r , shows that the 
second term on the right hand side of Eq. (3.24) does begin at order fi 3 . 

It now remains to use the explicit results for U r given in Eq. (3.23) and N r 
given in Eq. (3.14) and evaluate C for large fi keeping terms of order fi 3 and 
fi 2 . After some calculation it can be shown that the order fi 3 terms cancel 
leaving 



with 



cosh# fc n 2 kT ( 1 2 



Ei = 12 V(-l) fc - 

~ V ' ìsinh 2 ^ cu \ 



+ 



k=1 i^sinh 9 k uj \sinh# fc sinh Q k j J 

xcos^ttA;^ , (3.31) 



E 2 = y(-l) h \ — 2tt 2 A:T cosh^ ì g . n /^N 

\ sinh cj sinh 2 9 k J V a; / 



fc=i 



(3.32) 



* = 1 + E £^ «a. (3.33) 



~ /^wsinhé'fc V uo 
where 

- ^ (3.34) 

has been defined to save a bit of writing. The continuum approximation of 
Sec. 2 can be regained by dropping ali of the terms in £1,2,3 that have arisen 
from the de Haas-van Alphen part of the thermodynamic potential. This gives 
the familiar linear dependence on temperature [4]. 

The numerical calculations of [5] showed that as the temperature got suf- 
ficiently small there was a significant departure from the linear temperature 
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dependence in the sperine heat. Once again we will show that this follows from 
our analytical method, and that the origin of this behaviour is in the de Haas- 
van Alphen part of the thermodynamic potential. To do this we will evaluate 
the asymptotic expansion of the three sums defined in Eqs. (3.31-3.33) for 
large values of (3u. Leaving out the technical details of this for brevity, we find 
the approximate forms 



16tt 2 costi 2 (f (2/1-1)) 

S2 16. cosh 2 (^(2/2-1))' (3 ' 36) 
PU (3.37) 



4 costi 2 (f (2/2 - 1)) ' 

with Jx defined as in Eq. (3.19). The results in Eqs. (3.35-3.37) can be checked 
against a numerical evaluation of the sums defined by Eqs. (3.31-3.32) and 
found to be accurate for f3uj ~ 10 and p, not too dose to or 1. Once 
(5uj ~ 100 the results become very accurate even for p, dose to and 1. Thus 
for T < cu/100, the simple expressions in Eqs. (3.35-3.37) become reliable 
approximations for £1,2,3- 

If we use Eqs. (3.35-3.37) in the expression for the specific heat in Eq. (3.30) 
the result can be shown to vanish. The de Haas-van Alphen contribution to 
the specific heat cancels the continuum approximation to the leading order we 
are working to. The specific heat therefore vanishes as T — > faster than T. 
This is completely consistent with the numerical results found in [5]. Because 
the de Haas-van Alphen approximation, as well as the asymptotic evaluation 
of the sums leading to Eqs. (3.35-3.37) neglect terms that are exponentially 
suppressed, we suspect that the specific heat vanishes like e~ (XUJ l T for some 
Constant a as T -» 0, but we have not been able to establish this in any 
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simple way. Further support for this belief follows from the result we are 
able to establish for the 1- dimensionai gas in Sec. 4.1 below. A more refined 
estimate of the sums, as well as the de Haas-van Alphen contribution would 
reveal the exact nature of the T — ► limit. 



1.4 
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0.8- 



0.6 
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T = 0.2lo 



T = O.Ioj 




Fig. 2. (color online) This plot shows C/(6iV) 2 / 3 plotted over a range of N for 
four sample temperatures T = 0.1o»,T = 0.2u,T = 0.3uj and T = 0.5ùj. As the 
temperature increases the oscillation amplitude decreases and the curves approach 
the continuum limit of 1 in the scaled specific heat. As the temperature decreases 
there are significant deviations from the result found from using the continuum limit, 
and for very small temperatures the specific heat starts to become vanishingly small. 



For T > cj/100, but stili small, the results in Eqs. (3.35-3.37) start to become 
less reliable. As a check on our results against the numerical ones of [5] we 
plot the specific heat as found from Eq. (3.30) to demonstrate the de Haas- 
van Alphen oscillations. This is shown in Fig. 2. As the temperature increases 
we do find, as expected, that the contribution from £1,2,3 becomes smaller 
exponentially, and the linear behaviour with temperature is regained. 
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4 One and two dimensions 

We will examine the thermodynamics of trapped Fermi gases in one and two 
spatial dimensions, since these cases may be of relevance as limiting cases 
of the 3-dimensional gas in some situations. Because the methods used are 
similar to those described above in the 3-dimensional case we will be brief 
here and only exhibit key results that show a difference with results already 
obtained. 

4-1 One dimension 

With D = 1, Z{(5) in Eq. (2.15) has simple poles at (3 = (5 k with (3 k stili 
defined by Eq. (3.3). It is easy to show that the contribution from the poles 
with k ^ to 3(f?) is 

The contribution from the pole at /3 — was given in Eq. (2.18). After some 
calculation, making the approximation /i^$> T, it can be shown that 

1 ~ (-l) fc cos(27rfc/VccQ 

The calculation is very similar to the D = 3 case, so details will not be given 
here. The continuum approximation for Q, called Qq, was given in Eq. (2.24). 
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(4.1) 



(4.2) 



Using Eq. (2.24) and Eq. (4.2) in the general expression for the average particle 
number N in Eq. (2.27) results in 

- sinh(^0) • (4J) 

The first terni on the right hand side comes from Q Q , and the second terni 
from Q r , so that the continuum approximation would yield simply 

[i~uN (4.4) 

in this case. Again for large particle numbers we expect ^ > w to be a valid 
approximation. For T <C uj we can again find the asymptotic form of the sum 
in Eq. (4.3) and obtain 



N ~ 



1 1 , f /W2/f 



H h - tanh { - — — - 2 

2 2 } 4 U 



(4.5) 



and this proves to be very accurate for low values of T. 



The internai energy can be found to be 

U~£ + — T 2 + — + U r , (4.6) 
2u 6u 24 K J 



with 

27r , si . fu sin ( 27r /e u/u;) 7r coshé^ . _ , . . 1 . , „. 

U r ~— y -1 fe <^ ^ — v - 7 7 + sin (2-Kka uo) } 4.7 

/? É1 V ; \u> sinh0 fe /3wsinh 2 4 ^ w ; J 1 ; 

and defined by Eq. (3.34). We can calculate the specific heat using Eq. (3.24) 
with the result 

C ^{ 1 + E >- 12 |}' < 4 - 8 » 



26 



The sums Si, 2,3 are the same as those given in Eqs. (3.31-3.34). The result 
in Eq. (4.8) is very similar to that found in Eq. (3.30) in the 3-dimensional 
case apart from the overall factor that can be recognized as the continuum 
approximation for the specific heat of the 1-dimensional gas. We can therefore 
conclude immediately that as T — > 0, the specific heat vanishes faster than the 
linear temperature dependence deduced from the continuum approximation. 




Fig. 3. (color online) This plot shows the specific heat plotted over a range of 
temperature for the 1-dimensional gas. As the temperature increases the specific 
heat approach the continuum limit that behaves linearly on the temperature. As 
the temperature decreases the result deviates substantially from the linear result 
and the curve decays exponentially fast as found in Eq. (4.9). 

Because of the simplicity of the D = 1 case, we can shed some light on the 
behaviour of the specific heat as T — > 0. We note that the continuum approx- 
imation for the specific heat in Eq. (4.4) is also a solution when the de Haas- 
van Alphen contributions to N are included, since the sin term vanishes for 
integrai N. The sum £2 defined by Eq. (3.32) vanishes as well if fi = ujN, and 
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we cari use the approximation for Si given by Eq. (3.35), valici for T <C uj to 
obtain from Eq. (4.8) the simple result 

C ._^ exp (_£). (4.9) 

This demonstrates that the specific heat does vanish exponentially fast as 
T — > 0, not linearly. The role of the de Haas-van Alphen part of the ther- 
modynamic potential is again responsible for this behaviour. As a check ori 
this conclusion we have plotted the specific heat as a function of temperature 
in Fig. 3. The results can be seen to be consistent with our analytical result 
in Eq. (4.9) as T — > 0. For larger values of the temperature the specific heat 
approaches the linear temperature dependence predicted by the continuum ap- 
proximation. This can also be seen from the expression obtained in Eq. (4.8) 
since as /3u> becomes small, the sums £1,2,3 start to vanish as a consequence of 
the hyperbolic functions present in the expressions. 



4-2 Two dimensions 



Using D = 2 in Eq. (2.15) we find 



1 00 1 

X(E) = -—J^—{Ecos(2irkE/Lu) - — sìn^hE/u)} 



(4.10) 



This results in 



27T/3 fc tì ^ 



fi cos(27fkE/uj) 
uj sinh 9 k 

1 cosh 9k 



sinh 9 k sinh 6 k 



sm(2nkE/uj) \ , 



(4.11) 
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with 9k defined by Eq. (3.34) if we make the same approximations as in the 
D = 3 case. The (— l) k factor in the summand, present for D = 1, 3, is absent 
here, a result that is true for any even dimension. 



The average particle number can be shown to be 

,2 1 Jl 



N 



Li" 1 TX" 
J— I L /V 

2lo 2 12 Qf3 2 co 2 r ' 



with 



2tx ^ e /i sin(27rA;/i/a;) 7r cosh6> 



r ^ Iai 



+ — 



COS 



(2?rA;/i/a;)} 



^ L a; sinh 9 k ' sinh 2 ^ 

the contribution coming from the de Haas-van Alphen part of Q. 
In the continuum approximation (dropping the terni in N r ) we find 



fi 



LO 



~ (2iV) 1/2 , 



(4.12) 



(4.13) 



(4.14) 



to leading order for large N . However the de Haas-van Alphen part of the 
particle number can be shown to lead to the step-like behaviour we saw pre- 
viously in the 3-dimensional case. In the low temperature limit (T a;) an 
asymptotic analysis of Eq. (4.13) can be used to show that 



N ~ 



1 


>" 


2 3 


Li 


2 


.to _ 


2 


.CO 



+ 2 



tanh{^^-M) 
\ 2 \to [coìJ 



(4.15) 



Solving this expression for li gives rise to the approximate analytical solution 



fL 

LO 



(2N) 1/2 + 



(4.16) 



(Again we remind the reader that the square brackets in Eqs. (4.15) and 
(4.16) denote the greatest integer function.) We therefore conclude that the 
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step-like behaviour found earlier for the chemical potential when D = 3 is 
also found for the 2-dimensional gas. We have checked this result by solving 
Eq. (4.12) numerically and found that the approximation in Eq. (4.16) becomes 
increasingly accurate as T is reduced. Because the results resemble the similar 
behaviour exhibited in Fig. 1 we will not show them here. The jumps in the 
chemical potential occur, for large N, when N ~ |£ 2 for integrai i. (The exact 
result is t;£(£+ 1) analogously to the 3- dimensionai case studied in Ref. [5].) 

It is straightforward to obtain expressions for the internai energy and sperine 
heat. Once again the results are similar to those found in the 3-dimensional 
case, this time shown in Fig. 2, so need not be exhibited in detail. As T — > the 
result for the sperine heat vanishes much faster than the linear approximation 
C ~ 7r 2 T(2A^) 1 / 2 / '(Su) predicted by the continuum approximation. 

5 Anisotropie 3-dimensional potential 

5. 1 Thermodynamic potential 

As in the isotropie case, we will apply the analysis used by Sondheimer and 
Wilson [18] for the de Haas-van Alphen effect to obtain the thermodynamic 
potential. This allows a clear separation of the terms responsible for the con- 
tinuum approximation and those that lead to oscillatory behaviour in ther- 
modynamic expressions. The calculation of the thermodynamic potential will 
be used in the next two sections to obtain the chemical potential, particle 
number, and the sperine heat. 

Because we have already outlined the details of the Sondheimer- Wilson [18] 
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method above, we will begin with the main result that expresses the thermo- 
dynamic potential as (see (2.12)) 



where 3(E) is related to the inverse Laplace transform of the partition function 
Z(0) for the gas with /i = 0. (See (2.5) and (2.8).) 

To evaluate *ò{E) in Eq. (2.8) we require the partition function. For the case 
of a free Fermi gas confined by a simple harmonic oscillator potential we have 
the result in Eq. (2.15) with D = 3 taken there. The heart of the calculation 
consists of using Eq. (2.15) in Eq. (2.8), closing the contour in the left-hand 
side of the complex piane, and evaluating the result by residues as in the 
isotropie case considered in Sec. 3. 

As explained above, and as should be evident from Eq. (2.15), the poles of the 
integrand of Eq. (2.8) depend on the ratios of the three oscillator frequencies. 
There will be a pole of order 5 at (3 = as well as poles along the imaginary 
axis. It is easy to see from Eq. (2.15) that there are poles away from the origin 
at P — Pj where 



with j = 1,2,3 and kj = 0, ±1,±2, .... The ratio of oscillator frequencies 
determines the order of these poles, and the order affeets the residue. If none 
of the oscillator frequencies are rational multiples of each other, then the 
poles for /3i,2,3 will ali be distinct and therefore simple. However if any of 
the trap frequencies are rational multiples of any of the others, then some 
of the poles along the imaginary axis will no longer be simple. Therefore a 




(5.1) 
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complete analysis of the problem requires investigating a variety of different 
cases depending on the relative ratios of the three oscillator frequencies. In 
order to avoid too many different cases, and to stop the paper being overly 
lengthy, we will only look at the case of an axisymmetric trap characterized 
by 

uji = UJ2 = oj , CJ3 = Xoo , (5-3) 

where A is some real number that characterizes the anisotropy. There are then 
only two distinct cases to consider: A rational, and A irrational. The isotropie 
results can be recovered in the special case A = I. 

Regardless of whether A is rational or irrational, there is a pole of order 5 at 
j3 = for the integrand of Eq. (2.8). If we cali f2 the contribution to the 
thermodynamic potential that arises from this pole, then 



' r 3(8 + 3A 4 )^ 4 - 807T 2 TV(2 + A 2 ) + 1127r 4 T 4 ì , (5.4) 



5760Aa; 3 

as follows from the result in (2.26) with (5.3) used. The first term on the 
right-hand side of Eq. (5.4) can be recognized as the continuum (or semi- 
classical) approximation for the thermodynamic potential that arises from 
approximating the sums over the discrete energy levels with integrals, and is 
equivalent to results given in Ref . [4] . The other two terms are the next order 
corrections to the leading order behaviour. As we will discuss below, even in 
situations where the continuum approximation is valid, the leading order term 
may not be accurate if the trap becomes too anisotropie. If A becomes large, 
it is possible for the next to leading order terms, as counted in terms of /i, to 
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become non-negligible. 



5.1.1 A irrational 



When A is irrational it is easy to see from Eq. (2.15) that the integrand of 
Eq. (2.8) has doublé poles at 
2nik 



(3 = fa 



(5.5) 



and simple poles at 
2nik 



P = fa 



Xuj 



(5.6) 



where k = ±1, ±2, ... in both cases. (Because A is assumed irrational here, the 
fa poles are ali distinct from the fa poles.) It is a straightforward calculation 
to show that these poles make a contribution to Ò(E), that we will cali "òr (E), 
of 



3r(E) ="J^E U2 ■ ( V\\ { E * m ( 27rk -) 

47r 2 [ k 2 sm(nkX) L V co/ 



+ 



uj Xuj i / £J\ i 

— H cot(7rA;A) cos 27TA; — ^ 

7TK 2 J V UJ I > 



Xuj ~ cos - ttA;) 



E 



(5.7) 



8tt 2 ^ k 2 sm 2 (nk/X) ' 
This can be used in Eq. (5.1), and again making the approximation T <C /x 
we find 



1 00 1 
~j 77 E , 2 • , ,^ • htì { [ l + ^ cot(7rA;A) + 6 k coth fc 



x cos ( 27TÀ;- ) + 2nk- sin ( 2ttà;- 



cos 



(^-nk 



±y 

A(3 ^ksi^^k/X) sinh(0 fe /A) ' 



(5.8) 
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where 

2jT 2 k 



(5.9) 



At this stage we have not made any assumptions about the magnitude of (3uj 
in relation to [3 fi. A simple observation that can be made here is that for 
small values of fico the presence of the hyperbolic functions in Eq. (5.8) tends 
to suppress the contribution of Vt r to Q. Thus the oscillatory part of Q would 
be expected to become less important as the temperature is increased (say 
T ~ cu), and more important at lower temperatures (T < cu). For values of T 
larger than cu the leading part of Q would be expected to be well described by 
the continuum approximation. given by Qq. This will be verified in the next 
two sections. 



5.1.2 A rational 



For rational values of A we will define, without any loss of generality, 



A = - , (5.10) 



where p and q are two integers that are relatively prime. (That is, p and q 
have no non-trivial common factors.) The isotropie case considered above is 
a special case of this, corresponding to p — q — 1. The poles of the integrand 
of Eq. (2.8) stili occur at values of (3 equal to fa and fa given in Eq. (5.5) 
and Eq. (5.6), but this time some of the values of k result in some of the fa 
poles coinciding with those of fa. It is easy to show that Z(f3) has triple poles 
at 2niqk/uJ for k = ±1, ±2, . . .; doublé poles at 2mk/co for k = ±1, ±2, . . ., 
but k 7^ ±q, ±2q, . . simple poles at 2nqki / (pcu) for k = ±1,±2, ... but 
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k j£ ±p, ±2p, After some calculation, assuming that /j>Tas before, it 

cari be shown that 



1 oo 

^4e{[ 



p 2 + 2q 2 



(3 ^ \-2uj 2 pk sinh(q6 k ) 2Apq 2 k smh(q6 k ) 



ti 2 , 1_ 

'2(,)2nk. Vsinh ( r 



+ 



2[3 2 u 2 pk^smh(q6 k ) ' sinh 3 (g# fc ) 
cosh(g^ fe ) 1 



) 



cos(2irqk/j,/u — nkp) 



2f3pqk 2 sinh 2 (g# fc ) 47r 2 pg 2 À; 3 sinh(g^)- 

r 1 . g0fccosh(g0 fc )-i . -» 
27ipqujk 2 ^smh(qe k ) smh 2 (q6 k ) J V y P/ 
1 



1 oo 



4-7T/3 ^ /e 2 sin(7TÀA;) sinh 6 k 
2nkfj J . 



{ [1 + 7TÀ cot(7rAA;) + # fe coth0 fe ] 



x cos(27ikfj J /u) H sin(27rA;/i/cj)} 



i oo 



cos 



(5.11) 



4/3^1 A;sin 2 (7rA;/A)sinh(^/A) ' 

oo " 

Here ^ denotes that terms with k = q, 2q, 3q, . . . are omitted from the sum, 



k=i 

oo ' 



and ^2 denotes that the terms with k = p, 2p, 3p, . . . are omitted. In the 

k=i 

isotropie case, where p — q — 1, these last two restricted sums both vanish, 



and we recover the result for the isotropie case above. 



5.2 Particle number and chemical potential 



Having obtained the thermodynamic potential we can use it to calculate ali 
of the thermodynamic properties of interest. In this section we will study 
the chemical potential and its relation to the particle number. This can be 
obtained by solving (2.27) for fi. 

If we split f2 = iì + £l r as in (2.17) with fi given by Eq. (2.7) and £l r given 
by either Eq. (2.11) or Eq. (2.14), we can write N — No + N r in an obvious 
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way to find 



N = 




24Acj 3 



(2 + X 2 )lo 2 - 4tt 2 T 2 



(5.12) 



By temporarily neglecting the contribution N r to N coming from Q r , we obtain 



valid for large values of XN. For ÀiV > 1 we have /j, ^> lo. The result found in 
Eq. (5.13) contains the continuum approximation of Ref. [4] and the next to 
leading order correction. We will see that the next to leading order terni can 
have a significant effect for some traps. 

The approximation in Eq. (5.13) has neglected any effects that arise from the 
de Haas-van Alphen part of the thermodynamic potential. As in the isotropie 
case discussed in Sec. 3, the de Haas-van Alphen terms can make a signifi- 
cant alteration to the continuum result for the chemical potential for certain 
temperatures, and lead to the step-like behaviour discovered in the numerical 
result s of Ref. [5]. Because of the different expressions for Q found in the last 
section depending on the rational or irrational nature of A, we will look at two 
separate cases. 

5.2.1 X rational 

With A defined as in Eq. (2.13) we find, 



(6A7V )V3 + 1 ( 2 + A 2 



4tt 2 T 2 



) 



(6AA0~ 1/3 , 



(5.13) 




q/j 2 p 2 + 2q 2 



pio 2 sinh(g6> fc ) 12pq sinh(q9 k ) 

n 2 q / 1 2 
j _ i i 



) 



sm{2nqkji/LO — nkp) 



p(3 2 L0 2 ^smh(q8 k ) sinh 3 (g6> fc ) 
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27rg/i cosh(g^ fc ) -, 
+ ^^h^ 

+ 2^£ fcsin(7rAfc)smhg fc { + ^ COth ^ 

2irkii i 
cos(2ir k fi /lo) \ 

LO ' 
1 ^' Sili - 7T/fc 



E -27Z^^/7u • ( 5 - 14 ) 



2/3Awj^ sin 2 (7rA;/A) sinh(# fc /A) 



For fi ^> lo the leading contribution would be expected to be given simply 
by the first terni in the first sum on the right hand side of Eq. (5.14). It 
is the presence of the trigonometrie functions in Eq. (5.14) that result in 
the step-like behaviour of the chemical potential. Because of the presence of 
sinh(g# fc ) in the summand, where 9 k = 2ir 2 kT/LO, it can be seen that the 
sums converge very rapidly if T ~ lo. Fot T dose to lo very good analytical 
approximations can be achieved by merely taking the first terni of the sums. 
Because the temperature dependence of the sums is dictated by the magnitude 
of exp(— 27r 2 qT/Lo), the oscillations, although present, would be expected to 
be negligible and unobservable for T dose to lo. It is only once T is reduced 
below lo that the oscillatory, or step-like, behaviour will become observable. 



For values of T < lo, it is possible to rewrite the sums that occur in Eq. (5.14) 
in a way that is much more rapidly convergent. (As T becomes smaller than lo, 
the sums in Eq. (5.14) converge much more slowly rendering their evaluation 
more difficult. In addition, it is not sufficient to retain the first term, or even 
few terms, and achieve a reliable analytical approximation.) One way to obtain 
a reliable analytical approximation for T small in comparison with lo is to make 
use of the Poisson summation formula, familiar from the well-known Jacobi 
inversion formula for the #-function [21]. It can be shown that 
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^2, sin(27r/cm — irkp) 1 

^ sinh(27r 2 ÒA;) = Art> t{m 



1 , 

— (m-p/2) 



1 



sinh 
+ 



(rh-p/2) 



i, 



E 



2 ^ 



cosh(A;/ò) + cosh ( -(m — p/2) 



2tt6 

1 , (5.15) 



1 



where 



m = m — [m] , (5.16) 
p = pmod2. (5.17) 

Here [x] denotes the greatest integer that does not exceed x. The result in 
Eq. (5.15) is sufficient to evaluate the leading terni on the right hand side 
of Eq. (5.14). Other sums that occur in Eq. (5.14) may be evaluated in a 
similar manner. The advantage of this transformation is that for small values 
of b (where for Eq. (5.14) we have b oc T/uj) the slowly convergent sum on 
the left hand side of Eq. (5.15) is converted into a rapidly convergent sum on 
the right hand side. In fact, for b <C 1, we can safely neglect the sum that 
occurs on the right hand side of Eq. (5.15) and obtain a reliable analytical 
approximation from the first two terms on the right hand side, giving a very 
simple approximation. It is worth noting that the appearance of rh and p 
as defined in Eqs. (5.16) and (5.17) on the right hand side of Eq. (5.15) is 
necessary to preserve the periodic structure in m and p found in the originai 
sum on the left hand side; it is incorrect to simply take m and p on the right 
hand side of Eq. (5.15). An alternate way of obtaining the analytical terms on 
the right hand side of Eq. (5.15) makes use of the Euler-Maclaurin summation 
formula. This will not give the last term on the right hand side of Eq. (5.15) 
as it neglects exponentially small terms. 

If we take b <^ 1, and neglect the sum on the right hand side of Eq. (5.15), 
it is possible to approximate the result in an even simpler way. It is easy to 
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show that 



E 



oo 



sin(27rA;m — nkp) 



{1 -p + 2([m + p/2] -m)} . 



(5.18) 



sinh(27r 2 6A;) 



([x] is the greatest integer function mentioned above.) If we keep only the term 
of order fi 2 on the right hand side of Eq. (5.14), that would be expected to 
dominate the result for large values of fi ^> u>, the following simple approxi- 
mation cari be found: 



For the isotropie oscillator {p = q = 1) this reduces to the result in Eq. (3.20). 
When T < cu, but not necessarily such that T <C uj, it is more reliable to use 
the result of Eq. (5.15). 

When the temperature is low enough so that Eq. (5.19) is a reliable approx- 
imation, it is possible to solve for fi in terms of N and obtain the following 
analytical result: 



Because of the greatest integer function in Eq. (5.20) the result shows the 
step-like behaviour found in the numerical studies of Ref. [5]. In obtaining 
this result we have approximated N r with Eq. (5.19) and only included the 
leading order continuum approximation, the order fi 3 term, from Eq. (5.12). 
This means that the result would be expected to break down in situations 
where the neglected terms can become comparable in magnitude to those 
retained. One case of direct physical significance occurs when we allow the 
trap anisotropy parameter A to become large. In this case the subleading 




(5.19) 




(5.20) 
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Fig. 4. (color online) This figure shows the chemical potential in units of co 
plotted as a function of the particle number N for a range of temperatures 
T = O.lu, 0.05w, O.Olw with the anisotropy parameter A = 20. The curves with 
T = 0.05u; and O.Olw have been displaced vertically by 1.5 and 3 units respectively 
for clarity. The almost straight line shows the continuum approximation (6AÌV) 1 / 3 , 
and the step function over the T = O.lu; curve shows the analytical approximation 
[(6ÀÌV) 1 / 3 + ì(2 + A 2 )(6AiV)- 1 / 3 + \] found from Eq. (5.21) with p = 20, g = 1. 
(The temperature correction is negligible here and has been omitted.) 

part of Eq. (5.12) can lead to a non-trivial correction to the standard result of 
fi/u ~ (6AÌV) 1 / 3 as obtained in Eq. (5.13). It is not difficult to see that a better 
approximation for the chemical potential, that takes the de Haas-van Alphen 
part of Q into account, is 



a 1 
~ 

US q 



g(eAiV) 1 / 3 (l + l (^2 + A 



4vr 2 T 2 \ 
2 * (6AÌV)- 2 / 3 



0J* 



P 



+ 



P 



(5.21) 



in place of Eq. (5.20). (Recali that A = p/q here.) For small A and T uj 
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the result in Eq. (5.20) is a good approximation. For larger values of A it is 
necessary to use the more accurate estimation Eq. (5.21). 

In Fig. 4 we show the transition from smooth oscillations to step-like be- 
haviour of the chemical potential for a range of N and decreasing temper- 
atures. (Other ranges of N show a similar picture.) The trend towards the 
step-like behaviour as the temperature is reduced, predicted by the result in 
Eq. (5.21) is evident. We have chosen a large anisotropy parameter A = 20 to 
illustrate how the leading order continuum approximation becomes more in- 
accurate as the anisotropy increases. Had we plotted the improved expression 
of Eq. (5.13) it would have been centered on the chemical potential curves; 
it therefore provides a better estimate for the chemical potential, although 
stili missing the step-like features that are present. For smaller values of the 
anisotropy parameter, the result (6AiV) 1//3 becomes more accurate, and the 
plots look similar to those presented in the isotropie case described earlier. 



5.2.2 A irrational 

For irrational values of A we find 



1 00 1 

N r ~ 7—^ , . . — { [vrA cot(7rAA;) + 6 k coth# fc ] sin(27rfc///w) 

2pu^' 1 ksm(irXk) smhé> fe L 

cos(27vkfi/uj) \ 

ijj * 

1 ^ sin ( ^ - rrk 



2f3\u]^ 1 sm 2 (7rk/\)smh(6 k /\) ' ^ 5 ' 22 ^ 

This expression is simply the last two terms of Eq. (5.14) with the restrictions 
on the sums removed. For irrational values of A, the confluence of poles in 
Z(j3) that gave rise to the triple pole cannot occur here, and this is why the 
first term of Eq. (5.14) is not present. The leading order contribution to N r is 
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proportional to fi rather than fi 2 as we found in the case of rational A above. 
This leads to quant itati vely different behaviour for the chemical potential. 

It has proven very difficult to obtain reliable asymptotic approximations for 
the sums in Eq. (5.22), even for the leading order part of N r because of the 
presence of the sm(n\k) factor in the summand. The techniques that prove so 
useful in the rational case do not work here, as they ali involve the conversion 
of the sum into an ordinary integrai at some stage, and this requires the 
integrai to converge. The sin(7TÀ/c) factor prevents this here. As a result, our 
results are less complete than for the rational case examined earlier. 



Because of the lack of analytical expressions when A is irrational we will be 
brief here and only mention the main conclusions. For values of A ~ 1 (we 
used A = \/3) we found that although the oscillations in \i were present, 
they only became noticeable at very low values of the temperature, typically 
T = 10~ 2 u. This is expected as a consequence of the de Haas-van Alphen 
contribution to the particle number involving fi, rather than fi 2 as in the 
rational case. In addition the amplitude of the oscillations is much smaller 
than in the rational case. A conclusion is that the continuum approximation 
for fi as given by the first term on the right hand side of Eq. 5.13, is more 
reliable, even for temperatures as low as T = O.lu;, than in the rational case. 
For larger values of A (we used A = a/401 to compare with A = 20 in the 
rational case) similar conclusions were found provided that the extra term in 
the continuum approximation as given in Eq. 5.13 was included. 
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5.3 Specific heat 



The internai energy U is defined in terms of the thermodynamic potential il 
by (2.28). The specific heat C is then defined in terms of U by (2.36). With 
Q = Qq + fl r it can be readily seen that 



4 2 



8\lu 3 A8\lu 3 



(2 + Y)u 2 - l2n 2 T 2 (5.23) 



in agreement with what was found in (3.25). As in the previous section, the 
result for U r , coming from taking Q = Q r in Eq. (2.28), will be different 
depending upon whether A is rational or irrational. The results for U r will be 
given explicitly in Secs. 5.3.1 and 5.3.2 below. 

Although the internai energy U has a simple representation as U — U + U r , 
the same is not true for the specific heat because of N being held fixed in 
Eq. (2.36). It is easy to see that Eq. (2.36) can be expressed as in (3.24). 
The second term on the right hand side of (3.24) frustrates any attempt to 
express C simply as C = C + C r with C r the correction to the continuum 
approximation coming from the de Haas-van Alphen part of Vi. In order to see 
the effects that the de Haas-van Alphen contributions make to the specific heat 
we will initially ignore them, and work out just the continuum approximation. 

We will define Co to be the continuum approximation for the specific heat 
where we take U = U and iV = N in Eq. (3.24). Using Eqs. (5.23) and (3.3) 
we find (for T < uo) 
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The conclusion based on the continuum approximation is that at low tempera- 
tures the specific heat varies linearly with the temperature. However, as in the 
isotropie case, and the numerical studies of Ref. [5], this behaviour turns out 
to be significantly altered once the de Haas-van Alphen part of f2 is included. 



5.3.1 X rational 

If we concentrate on the leading order contribution to the specific heat for large 
/i, it would be expected, based on Eq. (5.24), that we would find C r oc fi 2 . We 
will therefore work consistently to order fi 2 when we evaluate C T . If we use fl r 
as given in Eq. (2.14) for f2 in the definition Eq. (2.28) we obtain 



3 



™" nfi 2 cos(2nkfi/u) 

^! ^ 2 sin(7rAA;)sinh^ 1 ' ; 

Terms of order fi or less have not been kept in this expression. (Recali that " 
on the second sum means that the terms with k = q,2q, . . . are omitted and 
that A = p/q here. 9 k was defined in Eq. (2.12).) 

'8U r 



When we compute 



dT 



, that is needed to find the complete expression 



fJ,,OJ 

for C, it can be seen that there will be a term of order fi 3 present rather than 
fi 2 as expected; however, it turns out that the term of order /x 3 coming from 
the first term on the right hand side of Eq. (3.24) cancels with an identical 
contribution coming from the second term, leaving C oc fi 2 as predicted. After 
some calculation, assuming that fi is large, it can be shown that 

^ 7T 2 T^ 2 r E 2 , , 



44 



with 



- [ cosh(gfl fc ) 7r 2 gfc / 1 2 \ 1 

1 ÉìlsinhVfe) ^ V sinh (^fc) sinh 3 (g# fc )/J 

x cos ( ^ nq ^ - Trkp^j , (5.27) 

V V- f 1 a cosh(gg fc ) | . ^Trgfc/x \ 

^[smh(g4) smlr(g# fc ) J \ w / 

£3 = 1 + -^ E • r? tt COS — ^-Trfcp . 5.29 

It is noteworthy that the dependence on the restricted summation in Eq. (5.25) 
that could have occurred in the end result has cancelled. The continuum ap- 
proximation of Eq. (5.24) is regained if we set £i,E 2 and S 3 to zero; this 
is expected because this is equivalent to dropping the contribution from the 
de Haas-van Alphen part of fi. 

For very low temperatures it is possible to obtain reliable analytic approxima- 
tions for £1,2,3- One way to do this is to use the Poisson summation formula as 
we did in Sec. 3.1 to approximate N r . If we assume that qT <C cu, it is possible 
to derive the following approximations: 



1 + Si ^ o) a 7 q q) xx , (5-30) 

4vr 2 gcosh 2 (f (m-\[qm\- |)) 



(3 2 u 2 (m - \[qm] - f q 



87rgcosh 2 (m - \[qm] - §)) ' 



(5.31) 



E „ ^ (5 32) 

3 4gcosh 2 (f (m-\[qm\-^) " 

As before, [x] is the greatest integer function and p = p mod2. We have 
not written out the terms in Eqs. (5.30-5.32) that behave like exp(— a/T) 
for Constant a as T — > 0. It is now possible to see that as T becomes much 
smaller than uj/q the specific heat behaves like Texp(— a/T), demonstrating 
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analytically that there is a deviation from the linear dependence predicted by 
the continuum approximation. As T increases, the terms in £1,2,3 in Eq. (5.26) 
begin to decrease in importance and the linear behaviour starts to become a 
good approximation. 
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Fig. 5. (color online) This figure shows the specific heat plotted as a function of the 
particle number iV over a range of N. The different curves show the results for five 
different temperatures. The anisotropy parameter A = 20. 

The behaviour of the specific heat with particle number and with temperature 
are shown in Figs. 5 and 6 respectively. In Fig. 5 the oscillations in the specific 
heat are clearly visible. The amplitude of the oscillations starts to decrease 
very rapidly as the temperature is increased. Once T becomes greater than 
about 0.5u they ali but disappear. As the temperature is reduced there is 
the competing effect that the specific heat starts to vanish exponentially fast 
as discussed above and illustrated in Fig. 5. Thus there is really a narrow 
range of temperature when both the specific heat and the amplitude of the 
oscillations are large enough to be potentially observable. In Fig. 5 it can be 
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Fig. 6. (color online) This figure shows the specific heat plotted as a function of the 

temperature for A = 20 and N = IO 5 both fixed. The straight line indicates the linear 

n 2 N 2 l 3 T 

temperature dependence given by the continuum approximation C = - — — ^ — . The 
other curve shows the effect if the de Haas-van Alphen contribution is included. 
noted that once T is about O.lu; there is a dramatic fall-off in the specific heat 

from the expected linear behaviour that is fully consistent with the analytical 
expressions obtained above. 



5.3.2 A irrational 

As in the calculation of the chemical potential, we were less successful at 
obtaining simple analytical expressions here than in the case of rational A. 
Using Eq. (5.8) for fì r in Eq. (2.28) we find 
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Fig. 7. (color online) This figure shows the specific heat plotted as a function of the 
anisotropy parameter with iV = IO 3 held fixed. A = p here with p an integer. The 
result is shown for three different temperatures T = 0.15ùj,0.3o; and 0.5u>. 
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This expression can be contrasted with that found in the rational case in 
Eq. (5.25). The de Haas-van Alphen contribution to the internai energy begins 
at order fi 2 rather than fi 3 . This has a direct consequence for the specific heat. 
If we use Eq. (3.24) to calculate the specific heat to leading order we find that 
to leading order in fj, the specific heat is given by the continuum expression 
in Eq. (5.24). For irrational values of A the de Haas-van Alphen contribution 
cancel out of the specific heat to leading order. There is a contribution from 
the de Haas-van Alphen part of the internai energy at sub-leading order /x; 
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however we will not give the expression here as it has only a negligitele effect. 
6 Discussion and conclusions 

The mairi conclusion of this paper is that the origin of the step-like features 
found in the numerical calculations of Schneider and Wallis [5] have the sanie 
origin as the periodicity found in the de Haas-van Alphen effect and that they 
can be approximated by analytical means. The general approach of Sond- 
heimer and Wilson [18], so useful in the analysis of the de Haas-van Alphen 
effect, can be used to great effect here to obtain analytical results for various 
thermodynamic quantities. In addition, it is possible to see that within this 
approach the continuum approximation used by Butts and Rokhsar [4] corre- 
sponds to the neglect of an infinite set of poles in the Laplace transform of the 
partition function. It is also possible to recover the continuum approximation 
using the methods described in the general framework of Ref. [22]. 

In Sec. 3 we obtained approximate analytical results for the 3-dimensional, 
isotropie harmonic oscillator potential. In the case T < w, a very simple re- 
sult (Eq. (3.21)) was found for the chemical potential that clearly exhibited 
the step-like features found in Ref. [5]. This was used to evaluate the sperine 
heat, and it was shown (in agreement with Ref. [5]) that there were significant 
deviations from the linear temperature dependence predicted by the contin- 
uum approximation. The specific heat was seen to exhibit a periodic structure 
in the particle number. 

In Sec. 4 we studied the trapped Fermi gas in 1 and 2 spatial dimensions. 
For the 1-dimensional gas we were able to show that as T — > the specific 
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heat vanished exponentially fast in the inverse temperature (Eq. (4.9)). The 
2-dimensional gas was similar in many ways to the 3-dimensional case of Sec. 3. 

The sanie methods that were applied in the isotropie case can be used to 
examine anisotropie potentials. The technical details are more involved due 
to the structure of the poles in the inverse Laplace transform of the partition 
function as described above. In addition to a periodicity in the particle number, 
there can also be a periodicity when the trapping potential is altered. This 
latter effect was not found for the isotropie case and it is worth commenting on 
why this occurs, in contrast with what might be expected from the de Haas- 
van Alphen effect where the thermodynamics shows a periodicity in the inverse 
magnetic field strength. The difference between the two cases is related to 
the relationship between the chemical potential and the particle number. In 
both cases (de Haas-van Alphen and trapped Fermi gas) the periodic structure 
results from a trigonometrie dependence on fi/cu. (For the de Haas-van Alphen 
effect cu is the cyclotron frequency associated with the magnetic field strength.) 
In the de Haas-van Alphen effect, when fi is solved for in terms of the particle 
number the leading order contribution turns out to be independent of cu. Thus 
the trigonometrie functions that involve fi/ cu exhibit the familiar de Haas- 
van Alphen oscillations as the magnetic field is varied. For the trapped Fermi 
gas in an isotropie potential we found \i oc cu, so that fi/ cu is independent of cu, 
although there is a dependence on the particle number. We claim that this is 
an artifact of the simplicity of the isotropie potential, and that the situation 
for anisotropie potentials is more interesting. 

The treatment presented here has only been performed for the free Fermi gas. 
As already mentioned in the introduction, the neglect of interactions for a sin- 
gle component gas is a good approximation [4,5]. In Ref. [6] a two component 
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model was studied that allowed for an interaction between the two spin compo- 
nents. A numerical study showed that as the interaction strength is increased 
the step-like behaviour, like that we have been discussing, is increasingly sup- 
pressed. For the single component gas, we expect that an analysis similar to 
that given by Luttinger [23] in the de Haas-van Alphen effect could be used 
to show that both the amplitude and period of the oscillations are not af- 
fected by interactions to leading order. (This conclusion does not hold for the 
non-oscillatory part of Q, that we have called Q .) It would be of interest to 
study the Luttinger analysis for a multi-component Fermi gas in more detail 
to make contact with the results of Ref. [6]. 
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